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Abstract. The effective potential approach for composite operators is generalized 
to non-zero temperatures in order to derive the equation of state for pure SU(3) 
Yang-Mills fields. In the absence of external sources, this is nothing but the vacuum 
energy density. The key element of this derivation is the introduction of a temperature 
dependence into the expression for the bag constant. The non-perturbative analytical 
equation of state for gluon matter does not depend on the coupling constant, but 
instead introduces a dependence on the mass gap. This is responsible for the large- 
scale structure of the QCD ground state. The important thermodynamic quantities, 
such as the pressure, energy and entropy densities, etc., have been calculated. We show 
explicitly that the pressure may vary continuously around T c — 266.5 MeV, whereas all 
other thermodynamic quantities undergo drastic changes at this point. The proposed 
analytical approach makes it possible to control for the first time the thermodynamics 
of gluon matter at low temperatures, below T c . We reproduce the properties of the so- 
called "fuzzy bag" -type models through the presence of the mass gap in our equation 
of state. An analytic calculation of the gluon condensate is obtained as a function of 
temperature. 
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1. Introduction 

The prediction of a possible existence of the quark-gluon plasma (QGP) created in the 
relativistic heavy ion collisions is one of the most interesting theoretical achievements 
of Quantum Chromodynamics (QCD) at non-zero temperatures and densities. A fairly 
full list of the relevant pioneering papers is given in [Tj EJ [3]) and the present status of 
the investigations of the properties of QCD matter is described in [HE]. 

The equation of state (EoS) for the QGP has been derived analytically up to the 
order g 6 ln(l/g 2 ) by using the perturbation theory (PT) expansion for the evaluation 
of the corresponding thermodynamic potential term by term ([61 [7] and references 
therein). However, the most characteristic feature of the thermal PT expansion is its 
non- analytical dependence on the coupling constant g 2 , which means that PT QCD is 
not applicable at finite temperatures. The problem is not the poor convergence of this 
series [TJ [9] but rather the fact that a radius of convergence cannot even be defined; 
any next calculated term can be bigger than the previous one. This is an in-principle 
problem which cannot be overcome. From the strictly mathematical point of view, four- 
dimensional QCD at non-zero temperatures effectively becomes a three-dimensional 
theory. At the same time, three-dimensional QCD has more severe infrared singularities 
[TO] and its coupling constant becomes dimensional. It is as a consequence of this that 
the dependence becomes non-analytical when using the dimensionless coupling constant 
g 2 . One also needs to introduce three different scales, T, gT and g 2 T, where T is the 
temperature, in order to try somehow to understand the dynamics of the QGP within 
the thermal PT QCD approach. 

At present, the only practical method to investigate the problem is lattice QCD 
at finite temperature and baryon density, which has recently shown rapid progress 
([U 13 [UJ [T2J [131 EH US] an d references therein). However, lattice QCD, being a very 
specific regular izat ion scheme, is primary aimed at obtaining well-defined corresponding 
expressions in order to get realistic numbers for physical quantities. One may therefore 
get numbers and curves without understanding what the physics is behind them. Such 
understanding can only come from the dynamical theory, which is continuous QCD. For 
example, any description of the QGP has to be formulated within the framework of a 
dynamical theory. The need for an analytical EoS remains, but, of course, it should 
be essentially non-perturbative (NP), approaching the so-called Stefan-Boltzmann (SB) 
limit only at very hight temperatures. Thus the approaches of analytic NP QCD and 
lattice QCD to finite-temperature QCD do not exclude each other, but, on the contrary, 
should be complementary. This is especially true at low temperatures where the thermal 
QCD lattice calculations suffer from big uncertainties (H El [HI [TJl [131 El US]- On the 
other hand, any analytic NP approach has to correctly reproduce thermal lattice QCD 
results at high temperatures (see papers cited above). There already exist interesting 
analytical approaches based on quasi-particle picture [TTj fT^ fT^ [20l l2T| l22| l23| l2^| [25] 
(and references therein) to analyze results of SU(3) lattice QCD calculations for the 
QGP EoS. 
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The main purpose of this paper is to derive the NP analytical EoS for the gluon 
matter (GM), i.e., a system consisting purely of Yang-Mills (YM) fields without quark 
degrees of freedom. The formalism we use to generalize it to non-zero temperatures is 
the effective potential approach for composite operators [26]. In the absence of external 
sources it is nothing but the vacuum energy density (VED). The approach is NP from 
the very beginning, since it deals with the expansion of the corresponding skeleton 
vacuum loop diagrams. The key element is the extension of our initial work [27] to 
non-zero temperatures. This makes it possible to introduce the temperature-dependent 
bag constant (pressure) as a function of the mass gap. It is this which is responsible for 
the large-scale structure of the QCD ground state. The confining dynamics in the GM 
will therefore be nontrivially taken into account directly through the mass gap and via 
the temperature-dependent bag constant itself, but other NP effects will be also present. 
Let us note that the temperature-dependent bag constant within the thermal PT QCD 
has been introduced into the Gibbs equilibrium criteria for a phase transition [28] (see 
also [16J). The effective potential approach has been already used in order to study the 
structure of QCD at very large baryon density for an arbitrary number of flavors [29] . 

The present paper is organized as follows. In section 2 the effective potential 
approach for composite operators is discussed in general terms. An explicit expression 
is obtained for the VED for pure YM fields. In section 3 the expression for the gluon 
pressure is derived taking into account the expression for the bag constant at zero 
temperature as a function of the mass gap [27]. Taken together this makes it possible to 
derive a formula for the pressure that is suitable for the generalization to non-zero 
temperatures. This consists of the two independent parts, describing the NP and 
PT contributions to the gluon pressure. In section 4 the generalization to non-zero 
temperatures is performed using the imaginary-time formalism. All the analytic results 
for the gluon pressure as a function of temperature are collected together in section 5. 
In section 6 the expressions are given for the main thermodynamic quantities, such as 
the entropy and energy densities, the heat capacity, etc., as functions of the pressure. 
In section 7 we discuss all the numerical results for the NP part of the gluon pressure. 
Section 8 concludes our discussion. Some explicit expressions for the summation of the 
thermal logarithms over the Matsubara frequencies are present in appendix A, while in 
appendix B a scale-setting scheme of our calculations is formulated. 

2. The Vacuum Energy Density 

The quantum part of the VED is determined by the effective potential approach for 
composite operators [26]. It is given in the form of the skeleton loops expansion, 
containing all the types of the QCD full propagators and vertices (for its pictorial 
representation, see [2TJ ) . Each vacuum skeleton loop is itself a sum of an infinite number 
of the corresponding PT vacuum loops (i.e., containing the point-like vertices and free 
propagators). Thus the effective potential approach makes it possible to calculate the 
VED from first principles, i.e., using only fundamental constituents of QCD: gluons and 
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quarks and their interactions or, more precisely, their propagators and vertices of their 
interactions. The number of the vacuum skeleton loops is equal to the power of the 
Planck constant, h, so that a series expansion for the effective potential is nothing other 
than the semiclassical WKB loop expansion, which, in general, is an asymptotic series. 
It has been widely used in quantum field theory [30]. In QCD the instantons have been 
discovered by using just this method [31] (and references therein). 

It is instructive to begin with some general expressions at zero temperature in this 
formalism. The gluon part of the VED to leading order (the so-called log-loop level ~ h, 
whose infinite series is shown in Fig. 1) is analytically given by the effective potential 
for composite operators [26] as follows: 

V(D) = l -J ^LSpiHD^D) - (D^D) + 1}, (1) 

where D(q) is the full gluon propagator and D (q) is its free counterpart. The traces 
over space-time and color group indices are assumed. It is clear from this that the 
effective potential is normalized to the free PT vacuum to be zero, i.e., V(D ) = 0. We 
note that the YM bag constant has been calculated to this order in |27J. 
The two-point Green's function, describing the full gluon propagator, is 

D^{q) = -i {T^(q)d(-q 2 ; + £V(<?)} ^ ( 2 ) 
where £ is the gauge-fixing parameter and 

T^(q) = g^u - = g^u ~ L^q). (3) 

Its free counterpart D = D® u (q) is obtaining by replacing the full gluon 
propagator's Lorentz structure d(—q 2 ;^) in Eq. (2) by unity, i.e., 

D%(q) = -i{T^(q)+ZL^(q)}± (4) 

In order to evaluate the effective potential (1) with Eq. (2), we use the well-known 
identity [32] 

SplniD^D) = 8 x Alndet^D^D) = 32 ln[(3/4)rf(-g 2 ; + (1/4)]. (5) 
Going over to four- dimensional Euclidean space in Eq. (1), one obtains (e 9 = V(D)) 




Figure 1. Infinite series for the gluon part of the VED to log- loop level. 
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where the constant a = (3/4) — 2 In 2 and an integration from zero to infinity over q 2 is 
assumed. The VED e g derived in Eq. (6) is already a colorless (color-singlet) quantity, 
since it has been summed over the color indices. It also does not depend explicitly on the 
unphysical (longitudinal) part of the full gluon propagator due to the product (D^D) 
which, in turn, comes from the above-mentioned normalization of the free PT vacuum 
to zero. Thus it is worth emphasizing that only the transversal ("physical") degrees of 
freedom only of the gauge bosons contribute to this equation. As a consequence, in the 
effective potential approach to leading order there is no need for the ghost degrees of 
freedom from the very beginning in order to cancel the longitudinal ("unphysical") 
component of the full gluon propagator. This role is played by the normalization 
condition. However, this does not work for the higher order vacuum skeleton loops. 
In this case and beyond the PT at any gauge (i.e., in the general case) the cancelation 
of unphysical gluon modes should proceed with the help of ghosts [27] . 

An overall numerical factor 1/2 has been introduced into Eq. (1) in order to make 
the gluon degrees of freedom equal to 32/2 = 16 = 8 x 2, where 8 color of gluons times 
2 helicity (transversal) degrees of freedom (see Eqs. (5) and (6)). 

In what follows the Lorentz structure d(q 2 ) = d(q 2 ; £) will be called the full effective 
charge ("running") or the gluon propagator invariant function (its form factor), for 
convenience. For the generalization of Eq. (6) to non-zero temperatures, the most 
important thing is to introduce correctly the above-mentioned bag constant. 

3. The gluon pressure at zero temperature 

The vacuum of QCD is a very complicated confining medium and its dynamical and 
topological complexity means that its structure can be organized at various levels: 
classical and quantum (see, for example [33], EU [35] and references therein). It is mainly 
NP by origin, character and magnitude, since the corresponding fine structure constant 
is large. However, the virtual gluon field configurations and excitations of the PT 
magnitude due to asymptotic freedom (AF) [36] are also present there. 

One of the main dynamical characteristics of the true QCD ground state is the 
bag constant B. Its name has come from the famous bag models for hadrons [37, 38J, 
but its present understanding (and thus modern definition) has nothing to do with the 
hadron properties. It is defined as the difference between the PT and the NP VEDs 
[391 STJl SU S2] . We can symbolically write B = VED PT — VED, where the VED is NP 
but "contaminated" by the PT contributions (i.e., this is a full VED like a full gluon 
propagator). Rewriting B as B = VED PT - VED = VED PT - [VED - VED PT + 
VED PT ] = VED PT — [VED TNP + VED PT ] = — VED TNP > 0, since the VED is always 
negative. The bag constant is nothing but the truly NP (TNP) VED, apart from the 
sign, by definition, and thus is free of the PT contributions ("contaminations"). In 
order to consider it also as a physical characteristic of the true QCD ground state, 
the bag constant correctly calculated should satisfy some other requirements such as 
colorlessness and dependence on the physical degrees of freedom, finiteness, gauge- 
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independence, no imaginary part (stable vacuum), etc. 

In our previous work [27], we have already derived an expression for the bag 
constant, which satisfies the above conditions: 



Bym = 16 



'eff 



d 4 q 



Mi + 3^V)]-^V) 



(7) 



P 9 = B YM - 16 



ln[l + 3d(q 2 )} - ^d(q 2 ) + a 



(8) 



(2vr) 4 

where symbolically shown q 2 ^ is the effective scale squared, separating the soft momenta 
from the hard ones in the integration over q 2 . This is a general expression for any TNP 
effective charge in order to calculate the bag constant from first principles. It is defined 
as the special function of the TNP effective charge integrated out over the NP region 
(soft momenta region, < q 2 < q 2 ff)- 

Adding the bag constant to the both sides of Eq. (6), and introducing the gluon 
pressure P g = e g + Bym, one obtains 

d 4 g 

(27T) 

The next step is to establish the relation between the full effective charge d(q 2 ) 
and its TNP counterpart d TNP (q 2 ). In our previous work jl3] it has been proven that 
the full effective charge depends explicitly and regularly on the scale A 2 (the so-called 
mass gap) responsible for the NP dynamics in QCD, i.e., d(q 2 ) = (q 2 ; A 2 ). The above- 
mentioned symbolic subtraction at the fundamental gluon level can then be defined as 
[27J S3]: d TNP (q 2 ; A 2 ) = d(q 2 ; A 2 ) - d(q 2 ; A 2 = 0) = d(q 2 ; A 2 ) - d PT (q 2 ). In this way 
the separation between the TNP effective charge and its PT counterpart becomes exact, 
but not unique. On how to make this separation exact and unique at the same time 
see subsection below. Evidently, both the TNP part and its PT counterpart are valid 
throughout the whole energy /momentum range, i.e, they are not asymptotics. Let us 
also emphasize the principle difference between d(q 2 ) and d TNP (q 2 ). The former is the 
NP quantity "contaminated" by PT contributions, while the latter one, being also NP, 
is free of them. 

However, this is not the whole story yet. Since the PT part d PT (q 2 ) contains the 
free gluon Lorentz structure d°(q 2 ) = 1, it should be extracted explicitly as follows, 
where we have dropped the explicit dependence on A 2 : 

d(q 2 ) = d TNP (q 2 ) + d PT (q 2 ) = d TNP (q 2 ) + 1 + d AF (q 2 ), (9) 

and d AF (q 2 ) describes the part responsible for AF in the PT gluon effective charge. This 
procedure is necessary in order to maintain the normalization of the free PT vacuum 
to zero. In this connection, let us stress that the second equality in this relation does 
not imply any violation of AF in the full gluon effective charge d(q 2 ), provided by the 
PT effective charge d PT (q 2 ) in the first equality of the same relation. Extracting d° — 1 
explicitly, we thereby subtract the divergent contribution associated with the constant a 
in Eq. (6), and thus the above-mentioned normalization condition will be automatically 
satisfied. We now replace all the effective charges as follows: d TNP (q 2 ) = ot TNP (q 2 ) and 
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d AF {q 2 



a 



AF 



(g 2 ). The explicit expression for the effective charge responsible for AF, 



that's a AF (q 2 ), will be given in part II of our work, since its explicit expression is not 
used here. 

Substituting the decomposition (9) into Eq. (8), lengthy algebra leads to 



Pg — PnP + PpT — B Y M + PyM + PpT, 

i.e., Pnp = P>ym + Pym- In Eq. (10) E>ym is given in Eq. (7), while 



YM 



and 



PT 



-16 



-16 



d 4 g 
(2*)* 



(2tt)' 



ln[l + \o? s NP { q 2 )\ ~ \c?. NP {f) 



ln[l + 



3a AF (q 



4 + 3a TNP {q 2 



r 



AF 



(10) 



(11) 



(12) 



Pym in Eq. (11) depends exclusively on the TNP effective charge. Together with the bag 
constant (7) it forms the NP part of the gluon pressure (10). Ppt in Eq. (12) contains 
the contribution which is mainly determined by the AF part of the PT effective charge, 
though the dependence on the TNP effective charge is also present (it is logarithmically 
suppressed in comparison to the pure AF term). If the interaction is switched off, 
i.e., putting formally a AF (q 2 ) 



a 



TNP 



(q 2 ) = 0, then P g = in accordance with the 
initial normalization of the free PT vacuum to zero. Evidently, just the right-hand-side 
of Eq. (10) is to be generalized to non-zero temperature, on account of the explicit 
expressions (7), (11) and (12). 



3.1. Confining effective charge 

The only problem remaining is the explicit expression for the TNP effective charge. 
Evidently, it has to be consistent with our initial work [27], where the bag constant was 
evaluated at zero temperature, i.e., 

«T"V) «rv) = ^ (is) 

where the superscript "INP" stands for the intrinsically NP effective charge. Here 
A 2 = A 2 JW is the Jaffe-Witten (JW) mass gap, mentioned above, which is responsible 
for the large-scale structure of the QCD vacuum, and thus for its INP dynamics [44]. 
Let us note that how the mass gap appears in QCD has been shown in our work |43j . 

A few additional remarks are in order. In (45] it has been shown that the TNP 
part of the full gluon propagator as a function of the mass gap also contains a term 
that is regular at origin. It is for this reason that it is not uniquely separated from 
the PT gluon propagator, which effective charge that is always regular at origin. We 
distinguish between the INP and the PT effective charges not only by the presence of 
the mass gap, but by the character of the IR singularities as well (45J H6]. So only after 
the replacement (13) the INP effective charge is uniquely and exactly separated from 
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its PT counterpart, and the obtained expression for the bag constant (7) becomes now 
free of all the types of the PT contributions ("contaminations"). 

In [46] we have shown that the so-called INP gluon propagator is, by construction, 
purely transversal in a gauge invariant way. It exactly converges to the gluon 
propagator, whose effective charge is Eq. (13), after the renormalization programme of 
the regularized mass gap is performed. This result has been obtained in the most general 
way, i.e., without making any truncations/approximations/assumptions or choosing a 
special gauge. Hence the expression (13) is not an ansatz, and thus it is mathematically 
well justified. However, it should also be physically well justified. The problem is 
that the gluon equation of motion is highly nonlinear, so the number of independent 
exact solutions is not fixed a priori (this number may be even bigger, depending on 
the different truncations/approximations/assumptions and the concrete gauge choice 
made). They should be considered on equal footing from the very beginning. In our 
previous work [45J, at least the two different general types of exact solutions for the 
full gluon propagator have been found as a function of the regularized mass gap: the 
first is smooth at small gluon momentum, allowing for the gluons to acquire an effective 
gluon mass (the so-called massive solution). The second is singular in the q 2 — > limit, 
so that the gluons always remain massless (the so-called nonlinear iteration solution). 
Only this solution has an affective charge (13) after the renormalization is completed. 
For the detailed arguments which makes its choice justified from the physical point of 
view as well (i.e., emphasizing its confining nature) see our paper [27]. In addition, let 
us note that the above-mentioned massive solution is not confining. 



4. Generalization to non-zero temperatures 

In the imaginary-time formalism [7J, H7] , all the four- dimensional integrals can be easily 
generalized to non-zero temperatures T according to the prescription (note that there 
is already Euclidean signature) 

dg 
(2vr) 

i.e., each integral over go of a loop momentum is to be replaced by the sum over 
the Matsubara frequencies labeled by n, which obviously assumes the replacement 
q — > u n = 2nnT for bosons (gluons). In frequency- momentum space the INP effective 
charges (13) becomes 

< NP {q 2 ) = af(q 2 ,a£) = ar>W) = (15) 



/(§W T E. g 2 = q 2 + % 2 = q 2 + ^ = ^ + ^ c n = 2nvrT, (14) 



n=— oo 



n 



and it is also convenient to introduce the following notations: 

cu = ^, (16) 

a AF (q 2 ) = a AF (q 2 , u 2 ) = a AF \u 2 , u 2 ) , where, evidently, in all the expressions here and 
below q 2 is the square of the three-dimensional loop momentum, in complete agreement 
with the relations (14). 



The Non-Perturbative Analytical Equation of State for the Gluon Matter: I 



9 



Introducing the temperature dependence into the right-hand-side of the relation 
(10), we finally obtain 



P g {T) = P NP (T) + P PT {T) = B YM {T) + P YM {T) + P PT (T). 



(17) 



4-1. Derivation of Bym{T) 

It is convenient to begin with Eq. (7) for the bag constant. In frequency-momentum 
space the temperature-dependent YM bag pressure becomes 



(27T)* 



Ml + 3a2" J V J t4!)]--a 



TNP ( r , 2 , , 2n 



.(11 



After the substitution of the expression from Eq. (15), one obtains 

i3„ +oo r + 3 a 2 



- j3 +oo / 

Bym(T) = 16 / (2^5 T i _E o p ( 



CJ 2 + W 2 / 4 w 2 + u; 2 



where we introduced the notations: 



q 2 + 3A 2 = ^Jcu 2 + m' 2 ff , m' eff = V3A. 



(19) 



(20) 



One of the attractive features of the confining effective charge (15) is that it allows 
for an exact summation over the Matsubara frequencies (see appendix A). Substituting 
all our results of the summation into Eq. (19), dropping the /3-independent terms [7j, 
and performing almost trivial integration over angular variables, one gets 



8 r u <>ff 
B YM (T) = — j / dco 

TV 2 JO 



3 A2 1 1 

-A 

4 



2/rMn 



1-e 



1 - e~^ u 



.(21) 



The uj e ff, which is the three-dimensional analog of q e ff in Eq. (7), is discussed in 
appendix B, where numerical values are given. 

It is convenient to present the integral (21) as a sum of several terms 



6 



B YM (T) = --A 2 B? M (T) - B&(T) - B^ M (T) 



16. 



7T 



7T 



>(2) 



>(3) 



(22) 



where the explicit expressions of all the terms are given as 



B?l{T) 
b ym( t ) 



dco— ; -. 

o e?" - 1 



J eff 



J eff 



duo uj 2 In ^1 — e 



du co 2 In ( 1 - e~^' 



Here and below iV 
mean number [7]). 



(23) 



1) 1 is the Bose-Einstein distribution (the so-called gluon 
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4-2. Derivation of Pym{T) 

In frequency-momentum space Eq. (11) becomes 

r d 3 a +OQ 
Pym(T) = -16 J T 



ln[l + ^ P (q>. 



-a T s NP [ 
4 



2 2\ 



thus determining the temperature-dependent YM part of the NP pressure of Eq. (17). 
After substituting of the expression (15), one obtains 



where we introduced the following notations: 



u> 



/ q 2 + -A 2 = 



v/3 



4w 2 + u^ 



1 , 



(25) 



(26) 



4 V " ""e/J> -wj 2 

Comparing Eqs. (19) and (25) one can write down the final result directly. For this 
purpose, in the system of Eqs. (22)- (23) one must change the overall sign, replace u>' by 
Gj and integrate from zero to infinity. Thus, one obtains 

where 



YM 



(T) 



4 2 i(^) - P&(T) 



(3) 



(27) 



Pym{T) 



p(2) 
YM 



JO 



do; ,, 

00 



a; 



7T 



-T 



1 6 
dw w 2 ln(l -e~ Pu ) , 

^°°da; a; 2 ln(l-e-^). 



(28) 
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The bag pressure (22) and the NP YM part (27) in GeV 4 units as functions of 
temperature T are shown in Fig. 2. The scale-setting scheme of all our numerical 
calculations in this paper is presented in appendix B. 

5. The gluon pressure at non-zero temperatures 

Summing up all the expressions and integrals (22)-(23) and (27)-(28), the gluon pressure 
(17) becomes 



where 



P g (T) = P NP (T) + P PT (T), (29) 



Pnp 

(T) = B YM (T) + P YM (T) 
6 1 6 

= -A 2 P 1 (T) + ^T[P 2 (T) + P 3 (T) - P 4 (T)}, (30) 



and 



while 



oo 



Pi(T) = / d^^-p, (31) 



ePt»-l : 



P 2 (T)= / doo oo 2 In (l - e-P" 

P 3 (T)= / dww 2 ln 
J o 



P 4 (T)= ^ doo io 2 In (1 - e- p " ) . (32) 

Let us recall once more that in all these integrals = T^ 1 , u) e ff along with the mass 
gap A 2 is fixed (appendix B), while u' and Q are given by the relations (20) and (26), 
respectively. In the formal PT A 2 = limit it follows that u = u' = u and the 
combination PziT) + P-3{T) — P^{T) becomes identically zero. Thus the NP part (30) 
of the gluon pressure (29) in this limit vanishes. 

In frequency- momentum space the PT part (12) of the gluon pressure (29) is 



Q roo +°° 

P PT {T) = -- 2 / duu 2 T 
n z Jo 



111 1 + 4 + 3a ™^ 2 ,^) J "4° W^>W) 



where the trivial integration over angular variables has already been carried out. 
Unfortunately, one cannot perform analytically (i.e., exactly) the summation over the 
Matsubara frequencies in this integral, using the AF expression for a AF {u 2 ,u 2 n ) [361116]. 
At this stage the PT part (33) remains undetermined and will be numerically evaluated 
elsewhere. Note that, since it becomes zero when the interaction is switched formally off 
(a AF (u 2 , oj 2 ) = 0), both terms (30) and (33) and hence the gluon pressure (29) satisfy 
the initial normalization condition. 
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Confining dynamics in Eq. (29) is implemented via the bag pressure (22) and the 
mass gap itself. However, other NP effects are also present via P Y m(T) and Ppt(T), 
given in Eqs. (27) and (33), respectively, though in Eq. (33) they are logarithmically 
suppressed. 



6. Main thermodynamic quantities 

Together with the pressure P{T), the main thermodynamic quantities are the entropy 
density s(T) and the energy density e(T). The general formulae which connect them 
are [7] 



s(T) = d -^l 

e(T) = T (^p) ~ PiT) = Ts(T) - P(T) (34) 

for pure YM fields, i.e., when the chemical potential vanishes. Let us note that in 
quantum statistics the pressure P{T) is up to a sign equal to the thermodynamic 
potential ft(T), i.e., P(T) = -tt(T) > [7j. 

Other thermodynamic quantities of interest are the heat capacity cy{T) and the 
velocity of sound squared c 2 (T), which are defined as follows: 

m de(T) r fds(T)\ 

CviT) = ^f- = T [-&r)> (35) 

and 

W = = w (36) 

de(T) c v (T) 

i.e., they are defined through the second derivative of the pressure. The so-called 
conformity 

C(T) = ^ (37) 

mimics the behavior of the speed of sound squared (36) but without involving such a 
differentiation. 

The thermodynamic quantity of a special interest is the thermal expectation value 
of the trace of the energy momentum tensor. This trace anomaly relation measures the 
deviation of the difference 

e(T) - 3P(T) (38) 

from zero at finite temperatures, in the high temperature limit it must vanish. As a 
consequence it is very sensitive to the NP contributions to the EoS. It also assists in the 
temperature dependence of the gluon condensate [H] (see [19]), namely 

< G 2 > T =< G 2 > -[e(T) - 3P(T)], (39) 




Figure 3. The NP pressure, the entropy and energy densities, the heat capacity as 
functions of T/T c . The NP pressure has a maximum at T c = 266.5 MeV. Here the 
horizontal dashed line is the general SB number (40), while the vertical dashed line at 
0.6T C separates the low-temperatures region from the transition region (0.6 — 1)T C . 



where < G 2 >o=< G 2 >t=o denotes the gluon condensate at zero temperature, whose 
numerical value is discussed in appendix B. 

As mentioned above, the main purpose of this paper is to investigate analytically 
and calculate the NP part (30) of Eq. (29), while the corresponding investigation of 
the PT part (33) will be the subject of subsequent work. This will make it possible to 
complete the investigation, derivation and final calculation of the full GM EoS. 



6.1. SB limit 

The high-temperature behavior of all the thermodynamic quantities is governed by the 
SB ideal gas limit, when the matter can be described in terms of non-interacting massless 
particles (gluons). In this limit these quantities satisfy special relations [7J 



3P SB (T) _ e SB (T) _ 3s SB (T) _ c v(SB) (T) _ 24 _ 2 



T 4 T 4 4T s 4T s 4{ V«5.26,T->oo,(40) 

and 

C SB (T) = c\ SB) {T) = I T^oo, (41) 

on account of the previous relations and their definitions in Eqs. (36) and (37). The 
trace anomaly relation (38) also satisfies the SB limit, namely 

e SB (T) - 3P SB (T) = 0, T — )• oo, (42) 



as it comes out from the relations (40). In what follows its right-hand-side will be called 
the general SB number. 
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Figure 4. The NP velocity of sound, conformity and the trace anomaly relation as 
functions of T/T c . 



0.115 



0.110 



0.105 



(G 2 ) 



Figure 5. The NP gluon condensate (39) in GeV 4 units as a function of T/T c . It 
shows little temperature dependence below T c and grows rapidly above T c . Here the 
solid line is its value at zero temperature < G 2 >o (see appendix B). 



7. Numerical results and discussion 



All our numerical results obtained for the thermodynamical quantities discussed in the 
previous section and calculated with the help of the NP pressure (30) are shown in Figs. 
3, 4 and 5. The NP gluon pressure has a maximum at some finite ("characteristic") 
temperature, T = T c = 266.5 MeV. This means that it may change continuously 
its regime in the near neighborhood of this point in order for its full counterpart 
to achieve the thermodynamic SB limit (40) at high temperatures. However, this is 
not possible for the other NP thermodynamic quantities, as is seen from the curves 
shown in Fig. 3 (none of the power-type fall off around T c can be continuously 
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transformed into the almost constant behavior at high temperatures). In order to 
achieve the corresponding thermodynamic SB limits at high temperatures their full 
counterparts would have to undergo drastic changes around this point. As we already 
know from the thermodynamics of SU(3) lattice QCD [19] (see [50l [5TJ [52]) the energy 
and entropy densities have discontinuities (jump or, equivalently, step discontinuities) 
at a temperature of about T c = 260 — 270 MeV. Our characteristic temperature 
T c = 266.5 MeV is surprisingly very close to this value. If we were not aware of the 
thermal lattice QCD results then we would be able to predict them. Since we are aware 
of them, these lattice results confirm our expectation of a sharp changes in the behavior 
of the entropy and energy densities in the region where the pressure is continuous. 

One of the interesting features of our calculations is seen in Fig. 3. At the 
characteristic temperature T c the NP pressure, the entropy and energy densities satisfy 
a SB-type relations, namely 

3P NP {T C ) _ e NP (T c ) _ 3gArp(T e ) _ 12 2 ^ 

7? ~~T^~ 4T C 3 "45^ [A6) 

where, obviously, the right-hand-side is half of the general SB number (40). In turn this 
yields 

C NP (T C ) = 1 e NP (T c ) - 3P NP (T C ) = 0, (44) 

so that the NP trace anomaly relation approaches zero from below in the T — > oo (/? — > 
0) limit (see Fig. 4), i.e, satisfying the SB relation (42) in this limit as well. 

Since the NP entropy and energy densities satisfy the SB-type relations (43), we 
expect for their full counterparts to have step discontinuities at T c , as discussed above. 
The NP heat capacity and the velocity of sound squared do not satisfy the SB-type 
relations (43), as can be seen in Figs. 3 and 4. We expect, therefore, for the full heat 
capacity to have an essential discontinuity at T c , since, in general, [cv(T c )] _1 = 0, and 
hence cj.(T c ) = due to the relation (36). They are defined through derivatives of the 
entropy and energy densities in Eqs. (35)- (36), i.e., they involve the second derivatives 
of the pressure. That is a reason why these thermodynamic quantities are too sensitive 
to the dynamical structure of the GM in the near neighborhood of T c . 

Moreover, due to the SB-type relations (43)- (44), and on account of the SB relations 
(40), for example it follows that 

~P SB (T)-2P NP (T)] T=Tc = 0, {^ T [Psb(T)-2P np (T)]} t=Tc = 0. (45) 

Combining these relations again with the SB- type relations (43)- (44) it is possible to 
derive some other exact relations between different thermodynamic quantities at T c . 
Their importance will be explicitly shown in part II of our investigation, where they will 
be extensively exploited. 

From Fig. 3 it clearly follows that the temperature range for the NP thermodynamic 
quantities can be divided into the three different intervals: 
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(i) . The low-temperatures interval up to its upper bound of about 0.6T C . It is 
defined by the first maximum of the heat capacity which appears in this region before 
temperature reaches T c . This region is obviously dominated by the NP contributions 
(30) to the gluon pressure (29). The behavior of all the NP thermodynamic quantities 
are shown in Figs. 5, 6 and 7. We do not expect any serious changes in the behavior 
of their full counterparts in this region. However, whatever changes may occur they 
will be under our control, since the NP part (30) is exactly calculated. The important 
observation is that all the NP thermodynamic quantities which have been properly 
scaled go down exponentially as temperature approaches zero. This is a general feature 
of the behavior of all the full thermodynamic quantities far below T c ; otherwise their 
zero temperature limit would not be realized. The almost constant behavior of the NP 
velocity of sound squared and conformity seen in Fig. 7 is due to the fact that they 
are the ratios of the corresponding thermodynamic quantities, see Eqs. (36)-(37). It 
seems that for the first time it is possible to predict the behavior of all the important 
thermodynamic quantities in this interval. We shall refer to the GM in this region as a 
"confining phase" (there are no lattice data for this region at all). 

(ii) . The transition interval (0.6 — 1)T C . In this region all the thermodynamic 
quantities should undergo sharp changes in their regimes, so that their full counterparts 
will be able to achieve their corresponding SB limits at high temperatures. The full GM 
pressure remains continuous in this region. However, the existence of maxima in the 
behavior of the NP thermodynamic quantities is expected in this interval, since they 
should "die" at high temperatures. The NP pressure is a part of the full GM pressure, 
which should be a continuously growing function of T throughout the whole temperature 
range. All other full thermodynamic quantities are continuous functions of T only below 
and above T c . Within our approach they are expected to have discontinuities of different 
character at T c , as mentioned above (apart from the gluon condensate). Following the 
authors of |53j, we shall refer to the GM in the transition region as a "mixed phase". 
The changes in the dynamical structure of the GM just in this region will determine the 
nature of the phase transition at T c . 

(Hi). The temperatures interval starting at T c . This interval itself can be clearly 
subdivided into two intervals, since in the integrals (31)-(32) u e ff = 3.75T C = 1 GeV: 

(a) The moderate-temperatures (1 — 3.75)T C , when the NP effects are still significant and 

(b) The high-temperatures starting at 3.75T C , when the NP effects become small. 
Therefore, even above T c , the GM can be understood as being in the two different 

forms. The first one can be considered as strongly coupled GM, where the NP effects 
are still important. The second one can be considered as weakly coupled GM, where 
the NP effects become already negligible. The typical temperature which is of about 
3.75T C = 1 GeV for the GM, is too high to be reached even at LHC. For QCD matter 
this temperature should be substantially decreased, and therefore be accessible at RHIC, 
and especially at LHC. The general feature for the behavior of all the properly scaled 
NP thermodynamic quantities above T c is their power-type fall off, while their full 
counterparts should show a power-type rise at high temperatures. Apparently, we may 
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refer to the GM in the moderate temperature region as an "extended mixed phase". 
The behavior of all the full thermodynamic quantities in both mixed phases is governed 
by the NP pressure (30) and by the first term in the PT pressure (33). We shall refer 
to the GM in the high temperatures region as an "AF phase", since the behavior of 
all the full thermodynamic quantities will be determined by the second term in the PT 
pressure (33) and the free gluons contribution. 

A few remarks in advance are here in order. From our discussion it clearly follows 
that the approximation of the gluon pressure (29) by its NP part (30) only (though 
exactly calculated in the present investigation) breaks down at T c , when the NP pressure 
achieves its maximum. For the derivatives of the pressure, this approximation breaks 
down even earlier, starting from about 0.6T C , as discussed above. This explains the 
sign change in the interaction measure and the supersonic value of the speed of sound 
squared both at T c as seen in Fig. 4. The rapid rise of conformity and the speed of 
sound squared above T c in Fig. 4 is also unphysical. This also explains why the gluon 
condensate in Fig. 5 increases drastically above T c . Taking into account the NP (30), PT 
(33), and free gluons contributions to the full GM pressure, all these unphysical effects 
will disappear in the full thermodynamic quantities, thus allowing them to achieve their 
corresponding SB limits. This will allow one to compare with lattice data above T c 
[HI |52] as well. 

7.1. The mass gap and a "fuzzy bag" term 

On account of the first of Eqs. (28) and Eq. (31), it is instructive to re-write the NP 
gluon pressure (30) equivalently as follows: 

P NP (T) = A 2 T 2 - - 2 A 2 P[(T) + - 2 T[P 2 {T) + P 3 (T) - P 4 (T)], (46) 
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Figure 7. The NP velocity of sound, conformity and the trace anomaly relation as 
functions of T/T c . They are shown in the low-temperatures region up to 0.6T C . 
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Figure 8. The NP trace anomaly relation scaled by T 2 in GeV 2 units (see Eq. (48)) 
as a function of T/T c . It clearly approaches a finite constant in the high temperature 
limit (compare with Fig. 4). 



where 

It scales as Tu e ff in the T — > oo limit, while the integrals P2CO, Ps{T) and Pa{T) 
remain unchanged, see Eqs. (32). Let us emphasize the explicit presence of the NP 
mass gap term in Eq. (46), namely A 2 T 2 . It will be present in the full EoS as well, 
since in the PT part of the gluon pressure the mass gap contribution is logarithmical 
suppressed (see Eq. (33)), and hence its contribution cannot be canceled. This is in full 
agreement with the so-called "fuzzy bag" model [52j [5U [551 [56] and with the massless 
Boltzmann stringy model [57J. The presence of the mass gap A 2 in our approach from 
the very beginning naturally explains the properties of such models. 
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Figure 9. The zero temperature gluon condensate and the NP trace anomaly relation 
with minus sign, as it enters Eq. (49). Both scaled by T 4 and shown as functions of 
T/T c . 




1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 

Figure 10. The NP gluon condensate scaled by T 4 as a function of T/T c . The zero 
temperature gluon condensate scaled by T A is also shown. At temperatures above 
1.5T C they go down as explicitly shown by Eq. (49). 



The NP trace anomaly relation scaled by T , i.e., 

*^M!xT' = fgg£3 - r f^' r > _> -const. T -» co(48) 

is shown in Fig. 8. It clearly goes to a finite constant (having the dimensions of A 2 ) in 
the high temperature limit, as it should in accordance with the previous discussion and 
the thermal lattice QCD results [I9| 1521 1551 158] . 
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This effect immediately shows up in the NP gluon condensate (39). Scaled by T 4 , 
analytically it becomes 

<G 2 > T <G 2 > [e NP (T)-3P N p(T)] 1 const. 

— — = — fg 7f 2 x 7jv ~ -^r^ T ~* °°-( 49 ) 

The difference in the fall off at high temperatures between the zero temperature gluon 
condensate and minus trace anomaly relation (as it enters Eq. (49)), both scaled by T 4 , 
is clearly seen in Fig. 9. Evidently, this difference is due to the asymptotic behavior of 
Eq. (48). At temperatures below T c there is no difference between < G 2 >t /T 4 and 
< G 2 >o /T 4 , which is as it should be (see, for example again Fig. 5). At temperatures 
above T c the difference between these two quantities is explicitly determined by Eq. (49) 
and is clearly shown in Fig. 10. The presence of the mass gap contribution in Eq. (46), 
and hence in the trace anomaly relation (38) scaled by T 4 , explains that the way it goes 
down is not 1/T 4 jlQl [58] but, as follows from our consideration, it decreases as 1/T 2 . 

The inclusion of the PT contribution [epriT)— 3Ppt{T)]/T a in Eq. (49) will produce 
a rather small numerical correction at high temperatures due to the SB limit (42) for 
the trace anomaly. The numerical results presented in Table 1 confirm this statement. 
It shows rather good numerical agreement between our NP analytical trace anomaly, 
as it enters Eq. (49), and its lattice counterpart at hight temperatures starting already 
from 2T C [32]. 

In conclusion, let us note that in lattice QCD at non-zero temperatures there is a 
possibility to distinguish between the magnetic (m) and electric (e) parts of the gluon 
condensate (see, for example [19], [591 EQ] and references therein). In [591 EQ] and [49J, 
such a gauge-invariant field-strength correlators in pure YM theory have been calculated 
at finite temperatures around T c and above T c , respectively. Following [60], we can 
formally write < G 2 >t=< G 2 >t + < G 2 m >t, but at this stage we do not know how 
to calculate these two parts separately within our approach. There is therefore as yet 
no possibility to compare with lattice data. In any case, such an analytical investigation 
at zero and non-zero temperatures is beyond the scope of the present investigation. 



Table 1. The NP analytical and lattice trace anomalies at high temperatures. 





Numbers 


T/T c 


Ours, Fig. 9 


Lattice, reference |49j 


2 


0.762 


0.860 


3 


0.453 


0.475 


4 


0.252 


0.255 


5 


0.166 


0.166 
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7.2. The NP dynamical structure of the GM. A brief description 

The NP dynamical structure of the GM within our approach is highly non-trivial. In 
the whole temperature range we have the two different massive gluonic excitations u' 
and u with the effective masses m! e ft = 1.17 GeV and m e ff = 0.585 GeV, respectively. 
Both effective masses are due to the mass gap A 2 , which is responsible for the large- 
scale structure of the QCD ground state. It is dynamically generated by the nonlinear 
interaction of massless gluon modes [4*3l H5j 06] . The first massive excitations can be 
interpreted as the glueballs, since Tn' e ^ is comparable to the masses of scalar glueballs 
(though it is lighter than the pseudoscalar ones) [61]. The second one rh e ff might 
be identified with an effective gluon mass of about (500 — 800) MeV, which arises in 
different approaches (see again the above-mentioned review (61] and references therein). 
We also have the two different massless gluonic excitations u: the NP massless gluons, 
propagating according to the integral P\iT) in Eq. (31), and almost SB massless gluons, 
propagating according to the integral P2(T) in Eqs. (32). We stress that the integral 
(31) should be multiplied by (6/7r 2 )A 2 , and all other integrals (32) are to be multiplied 
by (16/7r 2 )T, when one speaks about different NP contributions to the NP pressure (30). 

It is worth emphasizing once more that all these NP gluonic excitations are of 
dynamical origin, and only such can be accounted for within our approach. Whether 
the PT part (33) of the gluon pressure (29) will introduce a new massless or even massive 
gluonic excitations (or, equivalently, an effective gluonic degrees of freedom [HI [62]) is 
an open question at this stage. In our opinion, however, they may be really present 
in the full GM EoS. In addition to the mass gap the PT part (33) depends on the 
asymptotic scale parameter of QCD. The combination of these two fundamental scale 
parameters may lead to the creation of new massive gluonic excitations around T c , for 
example, pseudoscalar glueballs [60], EH [63] . In any case, we should definitely have the 
SB free massless gluons far away from T c . However, one thing should be made perfectly 
clear. Because of these and a possible new massive excitations in the GM, the behavior 
of all the thermodynamic variables in the moderate temperatures region (1 — 3.75)T C 
will be substantially different from the behavior of a gas of free massless gluons. This 
our qualitative yet conclusion is in agreement with recent lattice data and quasi-particle 
approaches both cited above. The recent lattice data will be reproduced quantitatively 
after the evaluation of the NP (30), PT (33), and free gluon contributions to the full 
GM pressure. 

The existence of the SB- type relations (43)- (44) indicates that in the transition 
region near T c a dramatic increase in the number of effective gluonic degrees of freedom 
will appear (for example, the glueballs will begin rapidly to dissolve). This effect will be 
reflected in the strong increase of the thermodynamic variables near T c , in agreement 
with recent lattice results [60, E3J EI]- This will lead to drastic changes in the structure 
of the GM. A change in this number is enough to generate pressure gradients, but not 
enough to affect the pressure itself. It varies slowly and therefore remains continuous in 
this region. At the same time, the pressure gradients, such as the energy and entropy 
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densities, etc., should undergo sharp changes in their behavior. The penetration of the 
PT part (33) of the gluon pressure below T c should not go deeply into the transition 
region (0.6 — 1)T C , since the INP effective charge is logarithmical suppressed. However, 
it should go sufficiently deep in order to eliminate "non-physical" maxima and other 
effects, created by the derivatives of Pnp{T) in the transition region. Instead of "non- 
physical" maxima, the full thermodynamic quantities (apart from the pressure) should 
have discontinuities at T c , as follows from our approach. This has also been established 
by the thermal QCD lattice cited above. 

Of course, not all the glueballs will be dissolved in the transition region. Some of 
them will remain above T c , along with other massive and massless gluonic excitations, 
forming thus a possible mixed and extended mixed phases around T c [53J. After the 
evaluation of the PT part (33), which is also "contaminated" by the NP contributions, 
we hope that the NP physics of the mixed phases will be well understood in the GM. At 
very high temperatures, starting at 3.75T C the NP effects become very small, and the 
structure of the GM will be mainly determined by the SB relations (40)- (42) between 
all the important thermodynamic quantities. 

8. Conclusions 

The effective potential approach for composite operators [26] has been generalized to 
non-zero temperatures in order to derive analytical EoS for pure SU(3) YM fields, shown 
in Eqs. (29)-(33). In its NP part (30) there is no dependence on the coupling constant, 
only the dependence on the mass gap, which is responsible for the large-scale structure of 
the QCD ground state. A key element of this work is the generalization of the expression 
for the bag constant at zero temperature [27] to non-zero temperatures. The NP part 
(30) of the gluon pressure (29) has been exactly evaluated, while its PT part (33) has 
been left undetermined at this stage, since it requires a separate investigation. 

Our main quantitative and qualitative results in this investigation, which will not 
be changed (or only slightly changed) after the inclusion of the PT part (33), are: 

(i) . In Eq. (29) the confining dynamics at non-zero temperatures (f5) is taken into 
account through the T-dependent bag constant (22) and the mass gap A 2 . 

(ii) . Other NP effects are also taken into account via the YM part (27). 

(iii) . The mass gap A 2 or, equivalently, u e ff is the only one independent input 
scale parameter needed to calculate the NP thermodynamic quantities. 

(iv) . The presence of the four different types of massive and massless gluonic 
excitations of the NP origin. 

(v) . The characteristic temperature T c = 266.5 MeV is a temperature at which the 
maximum of the NP part of the gluon pressure is achieved. 

(vi) . The low-temperatures region up to 0.6T C is under control. The exponential 
fall off in the T — y limit of all the main NP thermodynamic quantities (see Fig. 6) is 
expected to be preserved for their full counterparts as well. 
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(vii) . The behavior of all the thermodynamic quantities in the transition region 
(0.6 — 1)T C depends on how deeply the PT part (33) penetrates this region. 

(viii) . The existence of the SB-type relations (43)- (44) at T c shows that the 
structure of the GM below and above T c may be really rather different. 

(ix) . Since the NP entropy and energy densities satisfy them, we expect for their 
full counterparts to have jump discontinuities at T c , while the full pressure remains 
continuous. 

(x) . Since the NP heat capacity does not satisfy them, we expect for its full 
counterpart to have an essential discontinuity at T c , while the full speed of sound squared 
is to be zero at T c . 

(xi) . In the moderate temperatures region (1 — 3.75)T C the NP vacuum effects are 
still significant. 

(xii) . Because of this, our qualitative prediction is that in this region the behavior 
of all the thermodynamic variables will be substantially different from the behavior of 
a gas of free massless gluons. 

(xiii) . All the full thermodynamic quantities, therefore, approach their SB limits 
rather slowly. 

(xiv) . A possible understanding of new forms of the GM below and above 3.75T C 
as strong and weak coupling GM, respectively, has to be pointed out. 

(xv) . The existence of the mass gap term A 2 T 2 in the NP pressure (46), which 
remains in the gluon pressure (29) as well. 

(xvi) . Because of the mass gap term A 2 T 2 the full trace anomaly and the gluon 
condensate will go down as 1/T 2 at high temperatures, and not as 1/T 4 . 

In the subsequent paper (part II) we shall evaluate the PT part (33) of the gluon 
pressure (29), as well as include the free gluons contribution. It will make it possible 
to establish the order of the phase transition at T c . We will be able to compare our 
numerical results with thermal QCD lattice calculations (HI |52] at high temperatures 
above T c and in the transition region [52j |60j [63] as well. Only after completion of this 
programme, we will include the quark degrees of freedom in order to derive the NP 
QGP EoS within our formalism. For example, this will allow one to confirm a possible 
existence of Quarkyonic Matter (QM) [651 EE] and of a triple point in the QCD phase 
diagram (see [67] and references therein). Finally, let us note that recently [68J the 
pressure and the trace anomaly relation both scaled by T 4 have been calculated in 3d 
"electrostatic QCD" in apparent contradiction with the 4d lattice data [Hi [52]. It is 
worth noting once more that our calculations are in a good agreement with them (see 
Table 1 and discussion above). 
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Appendix A. The summation of the thermal logarithms 

In the second terms of Eqs. (19) and (25) the summation over the Matsubara frequencies 
can be done explicitly [7], as follows: 



E 



E 



n=-oo W + < n=-oo 



u 2 + (2nT) 2 n 2 \2n 



I) £ 



n 2 + (Pu/2tt) 2 



2tt ) f3u 



iPv-lJ 2u 



?/3u) _ 



l) ' 



In terms of the above- introduced parameters, the sums in Eq. (19) containing the 
corresponding logarithms look like: 

hi[3A 2 + to 2 + co 2 n } = lnu/ 2 + 2 ^ hi(2vr//3) > 2 + (/W /2tt) 2 ] (A.2) 



n=l 



and 
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n=— oo n=l 

It is convenient to introduce the notations: 
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In these expressions we introduced the following notations: 
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So the difference L(oj') — L(u) becomes 
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(A.7) 
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or, equivalently, 

iM -LM = --ln(^) + ln( — J. (A.8) 
From the relation (A. 6) it follows that 

so eq. (A.8) finally becomes 

L{u>) - L{u,) = - - In ^_ j + -W - u) + In [y— ^rj ■ (A.10) 

Appendix B. The scale-setting scheme 

Let us note that w e //, which appears first in the integral (21), is the only free parameter 
of our approach. In frequency-momentum space it is 



UJ eff = \]q 2 eff -uj 2 c1 (B.l) 

where we introduce the "constant" Matsubara frequency u c , which is always positive. 
Hence co e ff is always less or equal to the q e ff of four-dimensional QCD, i.e., u e ff < q e ff- 
One can then conclude that q e ff is a very good upper limit for possible values of w e //- 
In this connection, let us recall that the bag constant Bym at zero temperature has 
been successfully calculated at a scale q 2 e ^ = 1 GeV 2 , in fair agreement with other 
phenomenological quantities such as gluon condensate [27]. So u e ff is fixed as follows: 



^ff = = 1 GeV - (B-2) 

The mass gap squared A 2 , also calculated at this scale, has a value [27] 

A 2 = 0.4564 GeV 2 , A = 0.6756 GeV. (B.3) 

The effective gluon masses, defined in the relations (20) and (26), then become 

m' eff = v^A = 1.17 GeV, m eff = ^-A = 0.585 GeV. (B.4) 

The above-mentioned gluon condensate at zero temperature < G 2 > calculated at 
the scale (B.2) and at the same confining effective charge (13) in [27] is 

< G 2 > = (0\\g« u G%\0) = 0.1052 GeV 4 . (B.5) 

We need this value for the calculation of the temperature-dependent gluon condensate 
< G 2 >t via Eq. (39). Multiplying the right-hand-side of the relation (B.5) by (4a s /7r), 
where a s = 0.1187 |69| . it numerically becomes 



(Ol-^G^G^IO) w 0.016 GeV 4 . (B.6) 
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It is in a good agreement with its phenomenological value, < G 2 > p h= 
(0\{a s /7i)G a ^G a ^\0) « 0.012 GeV 4 , which can be changed by a factor of ~ 2, as 
mentioned in [3T] (see also |70j and references therein). In this connection, we 
recall [27] that the quark contribution to the bag constant, and hence to the gluon 
condensate, is approximately an order of magnitude less than the pure YM one. So the 
above-mentioned agreement of our YM value (B.6) with the phenomenological value is 
impressive. 
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